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ABSTRACT 

Recent studies have shown that at low metallicities Doppler-detected planet-hosting stars tend to have high a-content and to belong 
to the thick disk. We used the reconnaissance spectra of 87 Kepler planet candidates and data available from the HARPS planet 
search survey to explore this phenomenon. Using the traditional spectroscopic abundance analysis methods, we derived Ti, Ca, and Cr 
abundances for the Kepler stars. In the metallicity region -0.65 < [Fe/H] < -0.3 dex, the fraction of Ti-enhanced thick-disk HARPS 
planet harboring stars is 12.3±4. 1%, and for their thin-disk counterparts this fraction is 2.2+1.3%. Binomial statistics give a probability 
of 0.008 that this could have occurred by chance. Combining the two samples (HARPS and Kepler) reinforces the significance of this 
result (P ~ 99.97%). Since most of these stars harbor small sized or low-mass planets we can assume that, although terrestrial planets 
can be found in a low-iron regime, they are mostly enhanced by ff-elements. This implies that early formation of rocky planets could 
start in the Galactic thick disk, where the chemical conditions for their formation are more favorable. 
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1. Introduction 

Since the first discovery of an exoplanet around a solar-like star 
(lMavor&Ouelozlll995h . more than 750 exoplanets have been 
discovered and more than 230 planet candidates have been 
announced dBatalha et alj|20 12h . Despite the rather large num- 
ber, we still have to be content with a very small number of 
low-metallicity planet-hosting stars (PHSs) because of the well- 
established metal-rich nature of the PHSs, especially gian t PHSs 
(lGonzalezlll998l: Gonzalez et al ]l200lHSantos et al.ll200lL | 2004[ 



Fischer & Valentil|2005t iNeves et al.ll2009l Uohnson et al.l l201C : 
Sousa et al.ll201 lbt lAdibekvan et alj|2012bl) . Study of the chem- 



ical properties of metal-poor (iron-poor) PHSs is very important 
for understanding if there is any chemical composition require- 
ments for their formation. 

Recently, some studies have show n that iron-poor PHSs 
tend to be enhanced by g-elements (e.g. | Havwoodll2008l l2009t 
iKang et all l20lTl lAdibekvan et al.l 1201 2al) . This enhancement 
suggests that in a metal- poor regime most of the PHSs belong to 
the Ga lactic thick disk dHavwoo dl20"08l 120091: lAdibekvan et al] 
l2012ah . Ho wever, further interpreta tions of this result have been 
ambiguous. iHavwoodl (120081 120091) suggests that metallicity is 



* Table with chemical abundances is only available at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (138.79.128.5) or via 
http : //cdsarc . u-strasbg . fr/viz-bin/qcat? J/A+A/547/A36 



not the causative parameter that determines the presence of gi- 
ant planets around Sun-like stars, and there might be another pa- 
rameter linke d to galactocentri c radius (such as density of H2). 
Alternatively, lAdibekvan et al] (1201 2al) proposes that a certain 
chemical composition, and not the Galactic place birth of the 
stars, is the determinating factor for iron-poor PHSs to lie in th e 
high-a/thick-disk region (see also discussion in lGonzalezl2009l) . 

The iron content is used as a proxy for the overall metal- 
licity of planet hosts in most studies. However, iron is not 
the only abundant refractory element in the solar system. 
There are other fairly abundant elements (e.g. Mg and Si) 
with condensation temperatures comparable to iron (lLoddersI 
120031) that are very important contribut ors to the compo sition 
of dust in planet-forming regions (e.g. iGonzalezI 120091) . The 
c-enhancement of iron-poor PHSs is not unexpected, since 
a high [X/Fe] ratio means higher "total metallicity" and is 
moreover supported b y the theoretical studies using the core - 
accretion mo del (e.g. Ilda & Linl 12 004: Mo rdasini et al] 120121) . 
Interestingly, lAdibekvan et al] (1201 2al) show that not only stars 
hosting giant planets, but also Neptune- and super-Earth class 
planet hosts are also enhanced by <r-elements in the metal- 
poor regime. This indicates that despite the observed correla- 
tion between stellar metallicity ( [Fe/HD and presence of low- 
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Fig. 1. Comparison of the stellar parameters from B 12 to those derived in this work by traditional spectroscopic analysis technique 
(black dots) and with the TMCalc code (red crosses). 



iBuchhave et al.ll20l2l:lAd ibekvan et alJl2012bh . there is a special 
chemical composition required for their formation. 

IBuchhave et alJ (120121 hereafter B12) report spectroscopic 
metallicities of 152 stars harbori ng 226 exoplanet ca ndidates 
discovered by the Kepler mission dBorucki et alj|2010h . includ- 
ing 175 objects that are smaller than 4 R @ . With precise stel- 
lar parameters they also provide the unprocessed extracted spec- 
tra (only the wavelength region rangs roughly 260 A) used in 
their analysis. In this paper we have used these spectra and data 
from I Adibekyan et al.l d2012bl) to explore the possible chemical 
anomalies (i.e., a-enhancement) described above. 



2. The sample and atmospheric parameters 

Recently, B12 used the reconnaissance spectra obtained by the 
Kepler Follow-up Observing Program to derive the stellar pa- 
rameters for 152 stars (Kepler Objects of Interest - KOI) hosting 
226 planet candidates. These stars were chosen out of about 560 
KOIs (having reconnaissance spectra) as not very cool, r e ff > 
4500 K, and not fast rotators v sin; < 40 km s _1 , which allowed 
them to get more robust results. The spectra used in B12 were 
obtained with high-resolution spectrographs on medium-class 
telescopes, and some highest priority planet candidates (those 
of nea rly Earth-sized) w ere observed with the HIRES spectro- 
graph dVogt et al.ll 19941) on the 10-m Keck I telescope at Mauna 
Kea, Hawaii. The spectra have low/modest signal-to-noise ratio 
(S/N; S/N per pixel > 15). We refer the reader to B12 for more 
details. 

Our main goal is to study the chemical abundances of indi- 
vidual elements of planet-hosting dwarf stars with standard spec- 
troscopic analysis techniques using equivalent widths (EW) of 
the spectral lines. This technique requires higher S/N and lower 
v sin; than used in B12, where the stellar parameters were de- 
termined by matching an observed spectrum to a library grid of 
synthetic model spectra. For this reason we first selected only 
stars with v sin; < 10 km s _1 and \ogg > 4.0 dex, which re- 
duced the sample to 117 KOIs. Because many of these objects 
have multiple observations with different instruments, we were 
not limited in S/N (there is no numerical criterion), but during 
determination of chemical abundances we selected only those 
with reliable EW measurements (see Sect. [3j. Our final sample, 
for which we were able to determine elemental abundances, con- 
sists of 87 PHSs with 146 planet candidates. 



2. 1 . Testing the stellar parameters 

As mentioned before, in B12 the stellar parameters (T e s, logg, 
[m/H], and v sin;') were derived by cross-correlating the ob- 
served spectra against synthetic model spectra. To test how these 
parameters match those obtained with traditional techniques us- 
ing EWs, we selected the stars with highest S/N observed with 
Keck I (41 stars). For these stars we applied two different tech- 
niques to derive the stellar atmospheric parameters. 

In the first approach the spectroscopic stellar parameters and 
metallicities were derived based on the EWs of the Fe i and Fe n 
wea k lines by imposing excitation and ionization equilibrium 
(e. d Sousaet al. lE()TTbK . The spectroscopic analysis was com- 
pleted assuming local thermodynam ic equilibrium (LTE) with a 
grid of Kurucz atmosphere models dKuruczll 19931) . and making 
use of a recen t version of the MOOGQ radiative transfer code 
dSnedenlfl973l) . The quality of these spectr a was high enoug h 
to measure the EW using the ARES0 code dSousa et al.ll2007l) . 
Unfortunately, in the used spectral window there are 25 Fe i and 
5 Fe ii lines at most, which is not always enough to get precise 
parameters. 

In the second approach, we used the TMCalc code 
dSousa et al. 2012b. Th is code uses a recent line-ratio calibra- 
tion ( Sousa et al. 20101) to estimate the r e ff. Then, it uses a new 
direct spectroscopic calibration based on weak Fe i lines, which 
are expected to be less dependent on surface gravity and mi- 
croturbulence to estimate [Fe/H]. The code was combined with 
ARES to estimate both the spectroscopic stellar T & s and [Fe/H]. 
With this code we determined the mentioned parameters for 28 
KOIs. 

Comparing the stellar parameters derived with the two meth- 
ods we found strong discrepancies for the majority of the stars. 
This is not surprising because the first method is designed to use 
a long list of Fe i and Fe n lines. However, we found good agree- 
ment between the [Fe/H£| and determined with the TMCalc 
code and B12, except for three stars for which the differences 
in r e ff are higher than 250 K. The comparison of these two pa- 
rameters (except the three mentioned stars) is shown in Fig. Q] 
As can be seen, the two results agree well, although it seems 



at 



The source code of MOOG2010 can be downloaded 
http : / /www . as . utexas . edu/~chr i s/moog . html 

2 The ARES code can be downloaded 
http : / /www . astro . up . pt/sousasag/ar es 

3 The metallicity in B 12 is denoted as [m/H] and represents a mix of 
metals assumed to be the same as the relative pattern of the abundances 
in the Sun. 
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Fig. 2. [X/Fe] abundance ratios against [Fe/H] for the stars with 
and without planets from HARPS and Kepler surveys. The sym- 
bols and error bars (right panel) indicate the average and stan- 
dard deviation, respectively, of each bin (0.05 dex). The red 
squares and blue crosses represent stars with planets detected by 
HARPS and Kepler, respectively. The black dots refer to stars 
without a planetary companion. 



there is a small sy stematic trend at higher metallicities (see also 
iTorres et alJl2012h . The standard deviations of the differences of 
Teg and [Fe/H] between the two estimations are cr(T e g) — 82 K 
and <x([Fe/H]) = 0.07 dex (25 stars). This independent analysis 
also shows that the TMCalc code works well and can be very 
useful for determining precise parameters for large amounts of 
data. 

Assuming that the most reliable parameters were obtained 
for the stars for which our two methods gave similar results 
(lATeffl < 100 K and |A[Fe/H]| < 0.1 dex), we compared our 
atmospheric parameters (obtained with the first method) with 
those obtained in B 12. The comparison for these 16 PHSs is pre- 
sented in Fig. 03 As can be seen, the two results agree well. The 
standard deviations of the differences of the stellar parameters 
derived in B12 and with our first method are cr(T e ff) = 100 K, 
cr([Fe/H]) = 0.09 dex, and cr(logg) = 0.16 dex (16 stars). Based 
on this consistency we use the atmospheric parameters from B 12 
in our analysis here. 

2.2. Microturbulence 

The traditional spectroscopic abundance analysis methods us- 
ing the MOOG code and a grid of Kurucz model atmospheres 



requires four stellar parameters: [Fe/H], Teg, logg, and micro- 
turbulence velocity (f t ). B12 provides only the first three atmo- 
spheric parameters. Earlier studies of FGK dwarfs have shown 
that £ depends on T e g and log g (e.g. iNisse r3 ll98lHReddvetalJ 
12001 lATlende" Prieto et al. 20041) . Using the spectroscopic £ t of 



1111 FGK d warf stars from the high-S/ N and high-resolution 
HARPS data (ISousa et al.ll2QQ8Ll201 lblah . w 

ing expression 



we derive the follow- 



9.185 
-0.316 xlogg , 



x 10" 3 r e// + 3.436 x lO -7 ?^ 



(1) 



expecting an rms error of 0.18 km s . Here the £ is in km s , 
and r e ff and logg are in their traditional units. The stars in the 
sample have T e g ranging from 4500 to 6500 K, log g ranging 
from 3 to 5 dex, and [Fe/H] ranging from -1 .4 to 0.5 dex. We use 
this expression to calculate the £ for our Kepler sample stars. 

3. Refractory elements in planet-hosting stars 

The reconnaissance spectra used in this work cover from about 
70 to 260 A around A * 5200 A (see B12 for more details). In 
this spectral range we have four Ti 1 and one Ti n line, five Cr 1 
line s, and two Ca 1 lin es at most (for more details about the lines 
see iNeves et al ] l2009h . For these lines the EWs were measured 
using a Gaussian fitting procedure within the IRAlfl splot task. 
After measuring the EWs of spectral lines for the three elements, 
the chemical abundances were derived using an LTE analysis 
with the Su n as reference point with the MOOG spectral syn- 
thesis code (ISnedenll 19731) a nd a grid of Ku rucz ATLAS9 plane- 
parallel model atmospheres (lKuruczll993l) . The reference abun- 
dances used in the analysis were taken from lAnders & Grevessd 
(119891) . 

We derived the average abundances of the stars from all the 
available spectra taking their quality into account. Additionally, 
we removed the lines that were very different from the aver- 
age value. We are inclined to trust the results obtained from the 
spectra with higher S/N. For some stars there was only one low- 
quality spectra available with unreliable abundance results (very 
different abundances from different lines). We excluded these 
stars from our sample. At the end of this procedure, we finished 
with 75 KOIs with Cr abundances, 75 KOIs with Ca, and 86 
KOIs with Ti abundances. The final sample consists of 87 KOIs 
with at least one abundance result. 

In Fig. [2] we show the [X/Fe] abundance trends relative to 
[Fe/H] for the total Kepler sample. For comparison, the stars 
hostin g planets and field dwarfs without planets from the sam- 
ple of lAdibekvan et all (1201 2al HARPS sample) are also de- 
picted. In the same plot the averages and standard deviations of 
[X/Fe] ratios for each [Fe/H] bin (0.05 dex) are also displayed 
(right panels). These figures obviously show that the Kepler and 
HARPS PHSs have the same behavior, although the scatter of 
Kepler stars is much higher. We note that the methodology used 
for deriving the chemical abundances is exactly the same for 
both samples {Kepler and HARPS), and the errors due to the 
applied methodology (the errors induced by uncertainties in the 
model atmosphere parameters) are the same. This means that 
the heavy scatter found for Kepler PHSs is not likely to be astro- 
physical, and it indicates the low precision of the EW measure- 



4 IRAF is distributed by National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation, USA. 
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Fig. 3. The thin/thick separation histogram after subtracting the 
separation curve based on HARPS data {top). [Ti/Fe] against 
[Fe/H] for the metal-poor stars (bottom). The black dashed line 
indicates the separation between the thin- and thick disks. The 
blue crosses represent Kepler planet-hosting candidates. The red 
squares refer to the Jovian hosts and green triangles refer to the 
stars exclusively hosting Neptunians and super-Earths detected 
by HARPS. The black dots refer to the stars without a planetary 
companion. 

ments. The average scatter for KOIs in the metallicity region - 
0. 1 to 0.5 dex (in this regime the astrophysical scatter due to the 
Galactic chemical evolution is low) is about 0.1 1, 0.10, and 0.10 
dex for Cr, Ca, and Ti, respectively. For comparison the same 
scatters for HARPS PHSs are 0.02, 0.05, and 0.04 dex for Cr, 
Ca, and Ti, respectively, and for non-host stars these values are 
0.03, 0.04, and 0.04 dex. 

For the Kepler stars the line-to-line scatter errors we esti- 
mated are about 0.09, 0.10, and 0.10 dex for [Cr/H], [Ca/H], and 
[Ti/H], respectively. Taking the errors induced by uncertainties 
into account in the model atmosphere parameters (0.03, 0.04, 
and 0.04 for [Cr/H], [Ca/H], and [Ti/H]) and the errors in metal- 
licities (maximum 0.08 dex), we can expect the error of about 
0.13, 0.13, and 0.12 dex for [Cr/Fe], [Ca/Fe], and [Ti/Fe] ratios, 
respectively. 

4. Discussion and conclusions 

As already noted, the main goal of this work is to study the be- 
havior of Kepler PHSs at low metallicities, i.e., how these planet 
hosts are enhanced by <r-elements (relative to iron). In this study 
we used Ti as a proxy of overall a-content. We chose not to 
combine Ca with Ti, because the separation between the thin and 
thick disks based on the [Ca/Fe] is not as clear as for [Ti/Fel (e.g . 
iBensbv et al.ll2003t INeves et d.ll2009t lAdibekvan et at]|2012al) . 
In Fig. [3] (bottom), we p lot the distribution of stars both from 
lAdibekvan et al.l (1201 2bl) and KOIs from this sample in the 
[Ti/Fe] vs. [Fe/H] space. As can be seen, the stars are clearly 
separated into two groups according to their Ti content. The top 



panel of Fig. [3]shows the separation histogram after subtracting 
the thin - thick disk separation curve. We note that this separation 
is based only on HA RPS sample stars. 

As already found lAdibekvan et al. 2012b]) in the metal-poor 
regime, the frequency of PHSs in the chemically defined thick 
disk is higher than in the thin disk (f or details of the chemical 
separation see lAdibekvan et a"Dl201 ll) . Eight stars out of 65 are 
hosting planets (12.3+4.1%) above the separation line, and only 
three stars out of 136 harbour planets (2.2+1.3%) in the thin disk 
(low Ti-content). Because we have both the number of stars host- 
ing planets and stars without detected planets from the same sur- 
vey (HARPS), we can calculate the probability that 8 (or more) 
out of 1 1 PHSs found in this metallicity range belong to the thick 
disk (Ti-enhanced). The binomial statistics give a probability of 
0.008 that this could have occurred by chance. 

In the metallicity region -0.65 < [Fe/H] < -0.3 dex, there 
are six Kepler PHSs, five of which lie above the [Ti/Fe] sep- 
aration line (see Fig. |3j. If we assume that the Kepler targets 
also have a similar distribution in the [Ti/Fe] vs. [Fe/H] spaced, 
as we observed for HARPS (32% stars lie above and 68% bel- 
low the separation line), we can see that the probability of find- 
ing five (or more) out of six KOIs to be Ti-enhanced is about 
1.6%. Moreover, if we combine both samples (13 Ti-enhanced 
PHSs out of 17) the binomial probability gives a value of 3 x 
10~ 4 . This statistically significant result (3<x) suggests that in the 
metal-poor regime PHSs tend to have high Ti-content. 

Some of the Kepler planetary candidates are expected to be 
false positives (FPs) that do not turn out to be transiting plan- 
ets. Two stars in our metal-poor Kepler sample host three plan- 
ets, one hosts four planets, one hosts five, and only two stars 
host only one planet. The probability that all our planet candi- 
dates that are in the multiple systems are FPs is very low (e.g. 
iLissauer et al. 1120121) . and the FP probabilities for the two single- 
planet candidates are 2% for KOI105.01 and 4% for KOI373.01 
(iMorton & Johnsonll201 ll) . The probability that one of the KOIs 
or both of them hosts FPs is thus about 6%. If even one of the 
Kepler planet candidates is an FP, then we have four Ti-enhanced 
PHS out of five, for which the binomial statistics give a probabil- 
ity of about 4% that this could have occurred by chance. Taking 
the 6% FP probability in our sample into account we can calcu- 
late a binomial probability of about 98% that metal-poor Kepler 
planet candidates are really (not randomly) Ti-enhanced. The 
same statistics give a probability of about 99.96% (~ 3cr) for 
the combined (Kepler and HARPS) sample. 

The enhancement of metal-poor PHSs by Ti, suggests 
that these stars a r e also enhanced by oth er a-elements (e.g. 
INeves et all 120091: lAdibekvan et al.l l2012ah . including Mg and 
Si, which may be very important contributors to the composi- 
tion of dust in planet-forming regions and re present the prin- 
cipal components of ro cky-type planets (e.g. iGonzalej 120091: 
lAdibekvan etai1l2012bl) . 

B12 suggested that small planets may be widespread in the 
disk of our Galaxy, with no special requirement of enhanced 
metallicity for their formation. Interestingly, all six KOIs harbor 
small-size planets with a radius smaller than 4 R®. We note that 
six out of eleven PHSs (in the low-metallicity region) from the 
HARPS sample exclusively host Neptunes or super-Earth-class 
planets (less than 30 M e ), and only one lies in the Ti -poor region. 
The binomial statistics give a probability of about 4 x 10~ 4 that 
ten (or more) out of twelve small/low-mass planet hosts (com- 
bined sample) could randomly be Ti-enhanced. Considering the 



5 This can be very rough assumption because of the different fields 
of Kepler and HARPS. 
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possibility that Kepler planet candidates are FPs, the binomial 
statistics give a probability of about 99.95% (~ 3<r) that Ti- 
enhancement of small/low-mass planet hosts is not accidental. 
Our results suggest that, although terrestrial planets can be found 
in a low-iron regime, they are mostly enhanced by other metals. 
Moreover, high-a content at low-metallicities also hints that they 
belong to the Galactic thick disk, which implies that early for- 
mation of rocky planets could have started in the Galactic thick 
disk, where the chemical conditions for their formation are more 
favorable. 
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